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Abstract 

In this paper we consider the resonant frequency of split-ring resonators and double-bar 
resonators used to create artificial magnetic response at terahertz and optical frequencies. It 
is known that geometrical scaling of the resonant frequency of split rings breaks down at high 
frequencies (in the visible) due to electromagnetic properties of metals at those frequencies. 
Here we will discuss this phenomenon in terms of equivalent inductance and capacitance 
of the ring, derive an approximate formula for the resonant frequency of split rings in the 
visible and show how the limiting value of the resonant frequency of extremely small split 
rings and double bars depends on particle shape. 
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1 Introduction 



Realization of low-loss metamaterials possessing properties of an artificial magnetic in the 
visible would allow realization of artificial materials with negative index of refraction and 
open a way for many novel applications. One of the possible routes is the use of nano- 
sized resonant metal particles of complex shapes. Various shapes have been proposed, and 
the operational principle is basically the same: The particle shape is chosen so that one 
of the lowest eigenmodes has a current distribution with a large magnetic moment. In 
particular, several variants of split rings [1-3] and double bars [4,5] have been manufactured 
and measured in frequency range of hundrends of teraherts and in the visible. Also, arrays of 
pillar pairs have been studied in [6]. However, it has been noticed that the geometrical scaling 
of the particle resonant frequency breaks down when the working frequency gets higher, so 
that the resonant frequency saturates at the level of several hundreds of terahertz [7, 8] . 
This effect is explained by the plasmonic behavior of metals at optical frequencies, usually 
in terms of kinetic energy of electrons carrying current and an additional "electron self- 
inductance" [7]. 

In this paper we theoretically study the phenomenon of geometrical scaling breakdown 
and saturation of the resonant frequency using the equivalent circuit model. We show that 
the full model of the phenomenon requires not only an additional inductance, but also an 
additional capacitance. The model correctly gives the ultimate maximum resonant frequency 
which is equal to the plasma frequency of metal. This model allows us to study the depen- 
dence of the achievable resonant frequency on the particle shape. 

2 Magnetic mode resonant frequency in the visible 

In contrast to the microwave regime, in the visible range metals do not behave as ideal 
conductors, instead, their permittivity can be modeled by the Drude model: 



where oo p is the plasma frequency. This means that the fields in the bulk of metal strips 
cannot be neglected or modeled as fields in a thin skin layer. Let us assume that the 
current density is approximately uniform over the cross section of the strip and express the 
displacement current density in the metal strip as 



Next, we can relate the total current passing the ring cross section with the voltage along 
the ring perimeter, because the current amplitude equals / = Jdioh (again assuming uniform 
current distribution over the ring cross section) and the voltage is V = El c s. Here the 
effective ring length can be approximated as l c s = (ir/2)l, where / is the physical length, 
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Figure 1: The geometry of split rings and double-bar particles. The particles are fabricated 
from a strip of thickness h. The total length of metal strips in each particle is denoted by 
I. Arrows show the current direction of the magnetic mode. Magnetic flux and electric field 
exist not only in air but also inside metal strips. 



because the current distribution along an open loop or along short bars is approximately 
sinusoidal. The result reads: 
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This relation can be interpreted in terms of additional inductance and capacitance of the 
ring due to magnetic and electric fields inside the metal ring: 
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The additional capacitance and inductance can be then identified from the above equation 
as 
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Formula (5) for was derived earlier in [7] in terms of kinetic energy of electrons in metal. 
The resonant frequency of the split ring is given by 

uj = (6) 

y/(L + £ a dd)(C + C a dd) 

which differs from the earlier model in [7] by the additional capacitance term. 
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Figure 2: Equivalent circuit of a split ring or a double-bar particle in the visible. L add and 
C add model additional stored energy inside the metal volume. 



The parameters of the ideally conducting particle L and C depend on the dimensions 
and the shape of the particle. Considering the split-ring geometry, we can write: 



L ~ /i - hi 



wh 



(7) 



where we have approximated the split capacitance as the parallel-plate capacitance (5 is the 
split width). For the inductance estimations, we can use the formula for the round-loop 
inductance 

8a 

L = fi a In — (8) 

and estimate the equivalent loop radius a pa l e g/(2ir) = 1/4 and the equivalent wire radius 
r ~ (h + w)/4. Finally, the resonant frequency of split rings reads 
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Metal losses can be modeled by replacing formula (2) by 
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where 7 is the damping coefficient. The result for the resonant frequency (6) changes to 
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where the equivalent resistance reads 
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Figure 3: Resonant frequency of a split ring particle [formula (9)] as a function of geometrical 
scaling factor. Particle parameters are the same as in paper [8]. 

A numerical example of the dependence of the resonant frequency on the geometrical 
scaling factor is given in Figure 3. The particle is a split ring made of gold, and its dimensions 
are the same as in paper [8]. The plasma frequency of gold is assumed to equal u p = 13.8 x 
10 15 1/s, and the loss factor 7 = 107 x 10 12 1/s [9]. Despite several model simplifications, the 
estimated resonant frequency (316.4 THz) corresponds surprisingly well with numerical and 
experimental results of [8] (close to 330 THz). It is interesting to note that metal losses have 
negligible effect on the resonant frequency in this case. This can be understood by comparing 
the values of wn-^add and R, which differ only by replacement of the resonant frequency u>o 
by the loss factor 7. Because cun 3> 7, the contribution of the term proportional to R in (11) 
can be neglected. However, losses degrade the resonator quality factor and one has to take 
them into account especially at very high frequencies. In the optical region, the level of losses 
depends not only on the material but also on the particle shape, because of Landau damping 
(scattering on particle boundaries). This means that the choice of the optimal shape should 
include consideration of losses in addition to the resonant frequency of the particle. 

3 Discussion and conclusions 

In order to increase the resonant frequency, the denominator in formula (9) should be made 
as small as possible. Clearly, this has a fundamental limit since even in the limit of zero 
dimensions the limiting value is not infinity but equals to the plasma frequency of the metal 
(it is stated also in paper [8] that the limit is determined by the plasma frequency). This is a 
very expected conclusion because even if in the limit of small size the loop closes completely 
and there is zero magnetic flux outside the particle, electrons inside metal still can oscillate, 
and the resonant frequency of these oscillations is the plasma frequency. Inspecting formula 
(9) we can see that the preferred design should minimize the strip cross section area wh. 
This, however, has its own limit determined by increased losses for smaller values of wh. 
In addition, our results reveal that the limiting resonant frequency for particles with all 




the dimensions going to zero proportionally to the same scaling factor depends also on the 
particle shape: 

lim{ Wo } = -^= (13) 



[see formula (9)]. The ultimate highest resonant frequency equals to the plasma frequency 
of metal: ^ 

max{cj } = , = ^ P (14) 
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and it can be in principle achieved if in diminishing the particle dimensions 
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This design recipe suggests that, when reducing the particle size, the loop or bar length 
should be reduced in larger proportion than the gap width. This has a clear meaning of 
reducing both geometrical inductance and capacitance, since the inductance is proportional 
to I, while the capacitance is inversely proportional to 5. 

Finally, in choosing the particle shape it is important to remember that the current 
induced by external magnetic field (in the LC-circuit model) is proportional to uj 2 {C + 
C a dd)/(1— uj 2 Uq). This means that comparing two particles with the same resonant frequency 
ujq, the induced current is stronger in the particle with the larger total capacitance. This 
consideration suggests to decrease capacitance in a smaller proportion than inductance, 
if the goal is to design an optical metamaterial, but this contradicts with the criterium of 
raising the limiting resonant frequency (15). Thus, the particle shape choice is a compromise 
between the increased resonant frequency and strong magnetic response. 
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